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ABSTRACT 

We present the results of a high-resolution imaging survey of 22 brown dwarfs and very low mass 
stars in the nearby (~145 pc) young (^^1-2 Myr) low-density star-forming region Taurus-Auriga. We 
obtained images with the Advanced Camera for Surveys/High Resolution Channel on HST through 
the FbhbW iV), F775W (i'), and F850LP (z') filters. This survey confirmed the binarity of MHO- 
Tau-8 and discovered a new candidate binary system, V410-Xray3, resulting in a binary fraction of 
9 ±5% at separations >4 AU. Both binary systems are tight (<10 AU) and they possess mass ratios of 
0.75 and 0.46, respectively. The binary frequency and separations are consistent with low-mass binary 
properties in the field, but the mass ratio of V410-Xray3 is among the lowest known. We find that the 
binary frequency is higher for very low mass stars and high-mass brown dwarfs than for lower-mass 
brown dwarfs, implying either a decline in frequency or a shift to smaller separations for the lowest 
mass binaries. Combining these results with multiplicity statistics for higher-mass Taurus members 
suggests a gradual decline in binary frequency and separation toward low masses. The implication is 
that the distinct binary properties of very low-mass systems are set during formation and that the 
formation process is similar to the process which creates higher-mass stellar binaries, but occurs on a 
smaller scale. We combine the survey detection limits with models for planetary-mass objects to show 
that there are no planets or very low-mass brown dwarfs with mass > 3Mj at projected separation 
>40 AU orbiting any of the Taurus members in our sample, implying that planetary-mass companions 
at wide separations are rare. Finally, based on fits to the optical and near-infrared spectral energy 
distributions, we identify several BDs with significant (>1 mag) V-band excesses. The excesses appear 
to be correlated with signatures of accretion, and if attributed to accretion luminosity, may imply mass 
accretion rates several orders of magnitude above those inferred from line- profile analyses. 
Subject headings: stars:binaries:visual — stars:low-mass, brown dwarfs — stars :pre-main sequence 



1. INTRODUCTION 

Brown dwarfs (BDs) are objects with masses between 
those of stars and planets, insufficient to maintain sta- 
ble fusion reactions in their cores. BDs comprise a sig- 
nificant fraction of the total (sub)stellar content of the 
galaxy and are among our nearest neighbors (Reid et al. 
2004); in the past decade, field surveys have discovered 
several hundreds of BDs in the solar neighborhood (e.g. 
Delfosse et al. 1997; Kirkpatrick et al. 1999; Leggett et 
al. 2000). Soon after BDs were discovered, it was found 
that many, like stars, are members of binary systems. 
However, the properties of binary systems near and be- 
low the substellar boundary (Mprimary < 0.2Mq) ap- 
pear to be fundamentally different from those of higher- 
mass stars (0.3 < M„ri marv < 1.0M(n). Multi pHcity sur- 
veys of field T dwar fs llBurg-asser et al.ll2003t). L d warfs 
jKoerner et al."1999t [Close et al. ''200g: 'Bouv et al. 2003': 
ICizis ct al. 2003), and late M dwarfs (Si cglcr et al. 2005) 
have found lower binary frequencies (~15% vs 40-55%) 
and smaller binary separations (<20 AU vs < 1000 AU) 
than for field stars (Duquennoy & Mayor 1991; Fischer 
& Marcy 1992; Halbwachs et al. 2003). 

These results demonstrate that field binary prop- 
erties depend on mass. Unfortunately, binary fre- 
quencies for field stars are only reported for broad 
mass ranges, so they do not place strong constraints 
on the functional form of this dependence. Various 
groups interpret the transition in binary properties as ei- 
ther a sharp break near the stellar/substellar boundary 



ijKroupa et alJl2003HClose et al.ll200^ or a smooth mass 
deoendence l)Luhmanll2004bl ^. Also, since field BD sys- 
tems are typically old and possess lower binding energies 
than stellar binaries of equal separation, the results could 
be biased by the dynamical disruption of wide, low-mass 
systems. To begin testing this possibility, we performed 
a small multiplicity survey of low-mass stars and BDs in 
the nearby young OB Association Upper Scorpius (Kraus 
et al. 2005) and T Association Taurus (this work). In 
Upper Sco, we found several young binary systems, but 
the binary frequency, separations, and mass ratios were 
consistent with the field and somewhat (~2(t) inconsis- 
tent with higher-mass members of Upper Sco (Kohler et 
al. 2000). This suggests that dynamical evolution after 
the T Tauri stage probably does not produce the unique 
binary parameters of BDs; instead, the implication is 
that the mechanism by which binaries form depends on 
mass. 

Recent efforts to model low mass binary formation have 
typically assumed that a cluster of 5-10 protostellar em- 
bryos form from a single fragmenting cloud core (e.g. 
Kroupa 1995; Sterzik & Durisen 1998; Kroupa & Bou- 
vier 2003; Kroupa et al. 2003; Delgado-Donate et al. 
2003; Hubber et al. 2005); these embryos would then 
undergo dynamical evolution to form single stars and sta- 
ble multiple systems. However, the frequency of multiple 
stellar systems {Mprimary > 0.3Mq) in the field (35-57%; 
Duquennoy & Mayor 1992; Reid & Gizis et al. 1997) and 
in young associations (50-80%; Kohler et al. 2000; White 
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et al. 2006) has been interpreted by Goodwin & Kroupa 
(2005) to mean that the coUapse and fragmentation of 
a cloud core produces only 2 or 3 stars. Larger systems 
would eject more single stars and tight binaries than are 
observed. Reipurth & Clarke (2001) have suggested that 
evolution to a dynamically stable state could occur early 
in the formation process; the ejected embryos would then 
cease accretion and become BDs. 

The ejection process would preferentially disrupt wide 
BD binary systems, causing the deficit of wide systems 
seen in low-mass field binaries. However, simulations 
by Bate et al. (2003) find that the corresponding bi- 
nary frequencies and separations (<5% and <10 AU) 
are too low to be consistent with the field. Also, some 
wide very-low-mass binaries have recently been discov- 
ered in the field (Gizis et al. 2000; Phan-Bao et al. 
2005; Billeres et al. 2005). These systems are very 
weakly bound and most likely would not survive the 
ejection process, though recent simulations by Bate et 
al. (2005) suggest that they could form via simultaneous 
ejection of previously-unbound objects. Finally, some 
models predict that ejection would alter other proper- 
ties of BDs (spatial and velocity dispersion, disk lifetime, 
and accretion frequency). The preponderance of obser- 
vations show that these properties arc similar in the stel- 
lar and su bstellar reg imes ifWhite & B asri 20 03: Luhman 
|2p04b; W hite fc Hil lenbrand 2004; Mohantv et al. 2005); 
the strong similarity between the two regimes suggests 
that brown dwarfs for in a manner similar to stars and 
thus that BD binaries form like stellar binaries, though 
possibly on a smaller scale. The implication is that all 
binaries share a common formation mechanism, fragmen- 
tation of a single collapsing cloud core, and that this 
mechanism features a mass dependence that remains un- 
explained by theoretical models. 

Multiplicity surveys of the field and of nearby stel- 
lar populations (e.g. Kohler et al. 2000; Luhman et 
al. 2005a; White et al. 2006) have placed some con- 
straints on the form of this mass dependence. In par- 
ticular. White et al. (2006) studied the mass depen- 
dence of multiplicity in a speckle interferometry survey 
of the nearby (^145 pc) young (~l-2 Myr) T Association 
Taurus-Auriga. This surveyincluded objects from 1.5 
Mq to the substellar boundary and found that the sepa- 
ration distribution and mass ratio distribution functions 
were mass dependent. Their results for the binary fre- 
quency were inconclusive, but suggested a possible slow 
decline with mass. Better statistics will be required for 
very low mass binaries in order to confirm this trend. 
In this paper, we present the results of an complemen- 
tary imaging multiplicity survey of very low mass stars 
and brown dwarfs to near the planetary mass regime in 
Taurus-Auriga. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample Selection 

The Taurus- Auriga association has been the target of 
many recent wide-field surveys to detect new low-mass 
members (e.g. Briceno et al. 1998; Martin et al. 2001; 
Briceno et al. 2002; Luhman et al. 2003a). These sur- 
veys identified candidate members based on their loca- 
tion on an optical or near-infrared color-magnitude dia- 
gram, and membership was then confirmed spectroscop- 
ically via the detection of lithium absorption, excess Ha 



emission, or low surface gravity, all of which are indi- 
cators of youth. We selected all 18 confirmed Taurus 
members with spectral type later than M5.5 discovered 
by these surveys. Three additional targets (V410-Xray3, 
V410-Anonl3, and GM Tau) are previously known Tau- 
rus members which have also been confirmed to possess 
spectral types in this range (Strom & Strom 1994; White 
& Basri 2003). Our final target (LH0419-H15) was chosen 
from a survey for Hyades members by Reid & Hawley; it 
is the only member of their survey with a spectroscopic 
detection of lithium, and they classify it as a likely Tau- 
rus member with spectral type M7 based on its apparent 
youth and distance. We list these 22 targets and their 
discovery or confirmation references in Table 1; this was 
a complete list of known Taurus members with spectral 
types later than M5.5 at the time the observations were 
proposed (January 2003). Since our targets include both 
very low mass stars and brown dwarfs, we hereafter refer 
to them as very low mass objects, or VLMOs. 

2.2. Observations 

Our images were obtained with the Advanced Cam- 
era for Surveys/High Resolution Camera on the Hubble 
Space Telescope, which has a field of view of 26x29"and 
distortion-corrected pixel size of 27 mas pix~^. In Table 
2, we summarize the exposure time and epoch of obser- 
vation for each target (Program ID: 9853). Observations 
were made between September 2003 and January 2004 
with the filters F555W (V), F775W (i'), and F850LP 
(z') at two dither positions near the center of the detec- 
tor and with two exposures per position. Total integra- 
tion times were 510, 300, and 200 seconds, respectively. 
The F555W exposure times for the brightest objects were 
reduced to 350 seconds to allow for additional short ex- 
posures in F775W and F850LP, which were close to the 
saturation limit in the full-length exposures and were sat- 
urated in the case of MHO-Tau-5 and MHO-Tau-8. Sat- 
uration was permitted in the full-length images to allow 
for comparable sensitivities to faint companions at wide 
separations. We chose the V band to maximize angular 
resolution (diffraction limit Ores.v — 58 mas) and the i' 
and z' bands to maximize sensitivity to very low mass 
companions. 

The raw images were calibrated and distortion- 
corrected by the CALACS pipeline during on-the-fly- 
reprocessing (Mack et al. 2003). Some cosmic rays re- 
mained, but their morphologies were substantially dif- 
ferent from stellar PSFs, so they were easily identified 
by visual inspection. 

2.3. Data Reduction 

Potential point sources were identified with the IRAF 
task DAOPHOT/DAOFIND, and we measured aperture 
photometry and point-spread function (PSF) fitting pho- 
tometry for all objects in each field using the DAOPHOT 
tasks PHOT and ALLSTAR. We describe these data re- 
duction procedures in more detail in our companion pa- 
per on VLMO multiplicity in Upper Sco (Kraus et al. 
2005) and the DAOPHOT package is described by Stet- 
son et al. (1987). We report aperture photometry for all 
isolated objects and PSF photometry for all close bina- 
ries. PSF magnitudes were corrected to match aperture 
magnitudes based on results for the twenty apparently 
isolated Taurus members in our sample (Section 3.1). A 
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TABLE 1 
Observations 



Target 


Uate 


F555W 


: 

Exposure Times (s) 

F775W 


F850LP 


— 

Discovery 


CFHT-Tau-1 


2904.3 


510 


300 


200 


Martin et al. (2001) 


CFHT-Tau-2 


2904.3 


510 


300 


200 


Martin et al. (2001) 


CFHT-Tau-3 


2901.5 


510 


300 


200 


Martin et al. (2001) 


CFHT-Tau-4 


2857.5 


510 


300 


200 


Martin ct al. (2001) 


KPNO-Tau-1 


2881.2 


510 


300 


200 


Briceno ct al. (2002) 


KPNO-Tau-2 


2903.3 


510 


300 


200 


Briceno ct al. (2002) 


KPNO-Tau-3 


2901.4 


510 


300 


200 


Briceno ct al. (2002) 


KPNO-Tau-4 


2904.4 


510 


300 


200 


Briceno ct al. (2002) 


KPNO-Tau-5 


2884.2 


350 


300 /30 


200/20 


Briceno ct al. (2002) 


KPNO-Taii-fi 


3054.5 


510 


300 


200 


RrirpTin ct al {90C\'?S 


KPNO-Tau-7 


2901.6 


510 


300 


200 


Briceno ct al. (2002) 


KPNO-Tau-8 


28-50.5 


350 


300/30 


200/20 


Briceno ct al. (2002) 


KPNO-Tau-9 


3025.8 


510 


300 


200 


Briceno et al. (2002) 


KPNO-Tau-12 


3028.7 


510 


300 


200 


Luhman et al. (2003a) 


KPNO-Tau-14 


2897.5 


350 


300/30 


200/20 


Luhman et al. (2003a) 


MHO-Tau-4 


2893.5 


350 


300/30 


200/20 


Briceno et al. (1998) 


MHO-Tau-5 


3023.0 


350 


300/30 


200/20 


Briceno et al. (1998) 


MHO-Tau-8 


3028.5 


350 


300/30 


200/20 


Briceno et al. (1998) 


LH 0419+15 


2901.0 


510 


300 


200 


Reid & Hawlcy (1999) 


V410-Xray3 


3029.4 


350 


300/30 


200/20 


Strom & Strom (1994) 


V410 AnonlS 


2903.4 


510 


300 


200 


Strom & Strom (1994) 


GM Tau 


2850.6 


350 


300/30 


200/20 


White & Basri (2002) 



^Observation Date: JD minus 2450000. 



preliminary PSF for each filter was constructed from the 
20 well-sampled VLMOs that appeared isolated under 
visual inspection. Since all of the targets were located 
near the center of the chip and have similar tempera- 
tures, image distortion and target color should not be 
important. There was some variation in the PSF FWHM 
from target to target (±5%), which we attribute to small 
orbit-to-orbit changes in focus. We investigated this by 
dividing our sample into two groups, based on whether 
the target PSF appeared to be narrower or wider than 
the average PSF, and constructing new average PSFs for 
each group. We then re-ran PSF photometry, but the 
modest decrease in the residuals did not reveal any com- 
panions which were not previously identified. Since it is 
not possible to determine which PSF is appropriate for 
blended binaries, we proceed using only the average PSF 
for the entire group. This choice could lead to systematic 
errors in the calculation of binary properties; we discuss 
these errors in more detail in Section 3.4. 

Transformations to ground-based magnitudes {V, 
SDSS i', and SDSS z') were described in detail in Kraus 
et al. (2005). We found constant corrections that do not 
depend significantly on temperature or surface gravity: 
TO555— my = — 0.16±0.03, m775— rrij/ = -|-0.07±0.03, and 
1^850 — = -|-0.03±0.03. The transformed magnitudes 
and uncertainties are listed in Table 2. The statistical 
uncertainties correspond to either the photon noise (for 
aperture photometry; typically <0.01 magnitudes) or the 
goodness of fit (for PSF-fitting photometry; 0.01-0.06 
magnitudes for objects with significant fits). System- 
atic uncertainties in the magnitude transformations and 
aperture corrections are ~0.03 magnitudes. The pho- 
tometry calculated from short exposures was consistent 
with that from long exposures, so we report only the long 
exposures. 



3. RESULTS 
3.L VLMO Binaries 

The present survey found no fully resolved low-mass bi- 
naries. However, as we discuss in Kraus et al. (2005), one 
limitation in ALLSTAR-based PSF photometry is that 
binaries with very close {^X/D) separations are often 
not identified, even when their combined PSF is elon- 
gated at a high confidence level. DAOFIND, the task 
which identifies potential objects in the images, identi- 
fies point sources based only on the presence of a dis- 
tinct peak. Thus, automated photometry will be biased 
against the detection of very close binaries. This limita- 
tion can be overcome for known or suspected binaries by 
manually adding a second point source in approximately 
the correct location and letting ALLSTAR recenter it 
to optimize the fit. We have already used this method 
to recover photometry for the candidate binary system 
USco-109 in Kraus et al. (2005). 

In Figures 1-3, we illustrate this technique with con- 
tour plots of the known VLMO binary MHO-Tau-8 (dis- 
covered by White et al. (2006) and independently redis- 
covered here), the new candidate VLMO binary V410- 
Xray3, and the apparently single VLMO MHO-Tau-4 in 
the F555W, F77.5W, and'F850LP filters. Since the long 
exposures in i' and z' for MHO-Tau-8 were saturated, 
we show only the short exposures. Neither of the two 
candidate systems is obviously resolved and neither was 
consistently reported as a double-source by DAOFIND, 
but the PSF for MHO-Tau-8 is slightly elongated along 
the y-axis and that of V410-Xray3 is slightly elongated 
along the x-axis relative to MHO-Tau-4. 

In Figure 1, we present contour plots of MHO-Tau-8 
in each filter. The first column shows the original im- 
age, the second column shows the result from fitting 
with one point source, and the last column shows the 
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TABLE 2 

Photometry of Very Low Mass Objects in Taurus 



Name 



CFHT-Tau-1 
CFHT-Tau-2 
CFHT-Tau-3 

CFHT-Tau-4 

KPNO-Tau-1 

KPNO-Tau-2 

KPNO-Tau-3 

KPNO-Tau-4 

KPNO-Tau-5 

KPNO-Tau-6 

KPNO-Tau-7 

KPNO-Tau-8 

KPNO-Tau-9 

KPNO-Tau-12 

KPNO-Tau-14 

MHO-Tau-4 

MHO-Tau-5 

MHO-Tau-8 

MHO-Tau-8 A 

MHO-Tau-8 B 

LH 0419+15 

V410 Xray-3 

V410 Xray-3 A 

V410 Xray-3 B 

V410 Anon-13 

CM Tau 



23.427±0.026 
22.520±0.014 
21.654±0.009 

21.556±0.008 
24.063±0.044 
21.376±0.007 
20.239±0.004 
24.722±0.071 
19.690±0.004 
22.292±0.012 
22.068±0.011 
19.261±0.003 
24.918±0.084 
23.228±0.023 
20.735±0.007 
18.678±0.002 
17.595±0.001 
17.951±0.002 
18.241±0.030 
19.525±0.048 
21.835±0.010 
18.319±0.002 



22.175±0.012 
17.577±0.001 



18.612±0.003 
18.112±0.002 
17.970±0.002 

16.920±0.001 
19.572±0.004 
17.670±0.002 
16.959±0.001 
20.072±0.006 
16.226±0.001 
19.097±0.003 
18. 358+0.002 
16.147+0.001 
20.035+0.006 
20.781+0.009 
16.297+0.001 
15.246+0.001 
14.489+0.001 
14.453+0.001 
14.911+0.033 
15.611+0.060 
18.028+0.002 
15.046+0.001 
15.046+0.014 
20.032+0.446 
17.665+0.002 
15.169+0.001 



16.612+0.002 
16.307+0.002 
16.119+0.001 

14.951+0.001 
17.596+0.003 
16.074+0.001 
15.484+0.001 
17.897+0.004 
14.706+0.001 
17..301+0.003 
16.661+0.002 
14.863+0.001 
18.056+0.004 
18.998+0.006 
14.502+0.001 
13.732+0.001 
13.114+0.001 
12.976+0.001 
13.556+0.021 
13.934+0.013 
16.488+0.002 
13.629+0.001 
13.812+0.011 
15.658+0.028 
15.782+0.001 
13.908+0.001 



13.740+0.024 
13.754+0.021 
13.724+0.023 

12.168+0.020 
15.101+0.038 
13.925+0.022 
13.323+0.019 
14.997+0.033 
12.640+0.020 
14.995+0.028 
14.521+0.030 
12.946+0.018 
15.497+0.042 
16.305+0.085 
11.907+0.019 
11.653+0.028 
11.070+0.023 
10.862+0.018 



14.364+0.029 
11.544+0.018 



12.954+0.019 
12.804+0.019 



12.537+0.023 
12.762+0.021 
12.861+0.023 
11.008+0.019 
14.247+0.037 
13.241+0.027 
12.501+0.021 
14.025+0.037 
11.918+0.022 
14.197+0.038 
13.828+0.026 
12.367+0.019 
14.660+0.039 
15.483+0.096 
10.805+0.021 
10.923+0.032 
10.390+0.029 
10.140+0.020 



13.549+0.027 
10.817+0.021 



11.659+0.020 
11.586+0.017 



11.849+0.015 
12.169+0.017 
12.367+0.023 

10.332+0.016 
13.772+0.034 
12.753+0.020 
12.079+0.019 
13.281+0.031 
11.5.36+0.016 
13.689+0.036 
13.272+0.032 
11.985+0.020 
14.185+0.053 
14.927+0.092 
10.269+0.018 
10.567+0.020 
10.063+0.020 
9.726+0.016 



13.079+0.035 
10.446+0.017 



10.958+0.015 
10.632+0.016 



^Uncertainties are statistical only; systematic uncertainties due to aperture corrections and conversion to standard systems are ~0.03 magnitudes. 
^Near-infrared photometry is taken from the Two Micron All Sky Survey (Cutri et al. 2003). We quote total system magnitudes for the two 
binary systems. 



result from fitting witii two point sources. Tlie maxi- 
mum and minimum pixel values are also given to allow 
quantitative comparison of the residuals to the original 
images. The common position angle of the residuals in 
the single-source fit seems to imply that this elongation 
is a real effect, and not simply noise. Normal levels of 
jitter were reported in the observation log, so this is also 
unlikely to be a systematic effect. Unfortunately, there 
were no other bright objects in the field to serve as PSF 
references. However, as we summarize in Table 3, the 
system properties in each filter are consistent with pre- 
vious observations and expected orbital motion (Section 
5.3). Since MHO-Tau-8 was independently identified as 
a binary system by White et al. (2006), we regard it as 
a confirmed discovery. 

In Figure 2, we present similar contour plots of V410- 
Xray3, showing both the long and short exposures in 
i' and z'. The residuals from the single-source fit are 
roughly aligned and the jitter levels were normal, as in 
the case of MHO-Tau-8, but the separations of the resid- 
uals are marginally lower, implying a smaller separation 
or larger fiux ratio. ALLSTAR was unable to fit two 
point sources in the V image and the fit for the secondary 
in the i' images was not statistically significant, but the 
fit for the z' images appears consistent and statistically 
significant in both the short and long exposures. The 
number statistics for the i' and z' images are similar and 
the z' filter has lower resolution, therefore the superior 
fit for the z' images suggests that a possible companion 
may be much cooler and redder. We summarize the sys- 
tem parameters as calculated from each image in Table 
3, though the results for the i' images should be used 



with caution. 

The significant residuals in the double-source fit and 
the scatter in system properties suggests that even if 
V410-Xray3 is a binary, the measured parameters are 
not very reliable. Since there are no background stars 
for comparison and no other observations to support its 
multiplicity, we suggest its classification as a candidate 
binary. FoUowup observation to confirm its existence and 
properties should be a priority, since if confirmed, its 
small separation and corresponding short orbital period 
(~50 years for a circular orbit) could allow for a dynam- 
ical mass determination in less than a decade. Our sub- 
sequent analysis will consider V410-Xray3 as a candidate 
binary, but since its binarity has not been confirmed or 
disproved, our discussion will refiect both possibilities^. 

In Figure 3, we illustrate the typical results for single 
stars with contour plots of MHO-Tau-4 and the residu- 
als after fitting with a single point source. We could not 
obtain a statistically significant fit for two point sources, 
suggesting that there are no binary companions at sepa- 
rations >4 AU with mass ratio near unity. 

One target VLMO, KPNO-Tau-14, was reported as a 
possible double-lined spectroscopic binary by Mohanty 
et al. (2005). We did not detect any PSF elongation for 
this target. Since the components of SB2s have similar 
brightness, the detection limits we find in Section 3.3 
imply its separation is less than our inner detection limit 
(~4 AU). Since it falls inside our survey limits, we will 

^ Observations conducted in February 2006 with the Keck-II 
telescope and Laser Guide Star Adaptive Optics have confirmed 
this candidate and will be reported in a future pubhcation. 
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Fig. 1.— Contour plots of MHO-Tau-8 for all three filters. 
The first column shows MHO-Tau-8, the second column shows the 
residuals from fitting it with one source, and the last column shows 
the residuals from fitting it with two sources. For residuals, con- 
tours are drawn at the 90%, 50%, and 10% levels of maximum 
(solid lines) and minimum (dashed lines). The peak pixel value 
in each original image is shown; the positive and negative peaks 
of the residuals are reported as a percentage of the original peak 
value. 



not consider it as a binary in our discussion. 

3.2. VLMOs and Background Stars 

In Figure 4, we present an i' versus i' — z' color magni- 
tude diagram of the 20 apparently single VLMO targets 
(filled circles), the components of the two candidate bi- 
nary systems (filled circles with error bars), and all other 
objects which were clearly resolved and detected at the 
5(7 level in both filters (open circles). The candidate sec- 
ondary component of V410-Xray3 is well off the right side 
of the graph; its extremely red color {i' — z' — 4.374) is 
probably due to the high uncertainty in its i' fit (Section 
3.1); since the secondary has an approximate spectral 
type of M7.7 (Section 4.3), we expect it to possess col- 
ors similar to those of CFHT-Tau-3: i' — z' — 2, and 
thus i' = 17.5. Also shown are the average main se- 
quence at the distance of Taurus (Hawlcv ct al. 20Q2), 1- 
a nd 2-Myr isochrones based on the evolutionary models 
of lBaraffe et al.l l|1998j) . and a reddening vector based on 
the extinction relations reported in Schlegel et al. (1998). 
The location of the isochrone is determined by convert- 
ing the predicted Iq magnitudes to i' magnitudes using 
the color (/p — i'=0.40) derived as in Section 2 and by as- 
suming the i' — z' colors found from the S PSS field main 
sequen ce for older, more massive dwarfs bv lHawlev et al.l 
()2002t) . As we show in Section 4.3, the assumption of 
dwarf i' — z' colors for these targets appears to be valid. 
The DUSTY models (Chabrier et al. 2000) should not 
be required for these VLMOs due to their high tempera- 
tures at this young age, though we will use these models 
to determine limits on faint companions (Section 4.2). 

In Figure 4, the binary components are located above 
the SDSS main sequence and well above the background 
population, so they are most likely association members; 
based on their close proximity and the low density of 
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pixels 

Fig. 2.— Contour plots of V410-Xray3 for all three filters. 
The first column shows V410-Xray3, the second column shows the 
residuals from fitting it with one source, and the last column shows 
the residuals from fitting it with two sources. For residuals, con- 
tours are drawn at the 90%, 50%, and 10% levels of maximum 
(solid lines) and minimum (dashed lines). The peak pixel value 
in each original image is shown; the positive and negative peaks 
of the residuals are reported as a percentage of the original peak 
value. 



association members we conclude that these are physi- 
cally associated companions. All other detected objects 
fall well below the empirical main sequence and are most 
likely background stars. One target object, KPNO-Tau- 
12, appears to be somewhat underluminous compared to 
the other Taurus VLMOs. The discovery survey (Luh- 
man et al. 2003a) also found KPNO-Tau-12 to be un- 
derluminous for its assumed age and distance, so this is 
most likely a genuine feature of the system. 

3.3. Sensitivity Limits 

We determined detection limits as a function of dis- 
tance from the primary stars via a Mo nte Ca rlo simula- 
tion similar to that of Metchev ct al] l)2003f) . We used 
the IRAF task DAOPHOT/ADDSTAR and our average 
PSF to add artificial stars at a range of radial separa- 
tions and magnitudes to the fields of GM Tau, KPNO- 
Tau-2, KPNO-Tau-7, and KPNO-Tau-4, which represent 
the full range of brightness in our sample. We then at- 
tempted to identify the artificial stars with the ALL- 
STAR PSF-fitting photometry package. Since the one 
confirmed VLMO binary in our sample (MHO-Tau-S) 
has already been resolved in other observations, our goal 
for these simulations was to determine the performance 
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TABLE 3 
Binary Properties 



larget 


Date 


beparatioii 
(mas) 


Position 
Angle(deg) 


Filter 


Flux Ratio 
Am (mag) 




q 


MHO-Tau-8 




37 


311 


K 


0.02 








3029 


54 


269 


V 


1.28 








3029 


39 


272 


i' 


0.70 








3029 


40 


268 


z' 


0.378 








<3029>= 


44±8 


270±2 






0.6±0.1 


0.75±0.07 


V410-Xray3 


3029 


56 


192 


i'(long) 
z' 


5.0 








3029 


45 


205 


1.85 








3029 


46 


234 


i'(short) 


3.7 








3029 


43 


203 


z' 


1.55 








<3029> " 


44±2 


204±2 






1.7±0.2 


0.47±0.07 



^Observation Date: JD minus 2450000. 
''Wliite et al. (2006) 

System values for MHO-Tau-8 are determined from the mean and standard deviation for all HST observations; values for V410-Xray3 are 
determined from the long and short z' observations 
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Fig. 3.— Contour plots of MHO-Tau-4 for all three filters. The 
first column shows MHO-Tau-4 and the second column shows the 
residuals from fitting it with one source. For residuals, contours are 
drawn at the 90%, 50%, and 10% levels of maximum (solid lines) 
and minimum (dashed lines) . The peak pixel value in e3x:h original 
image is shown; the positive and negative peaJcs of the residuals 
are reported as a percentage of the original peak value. 



in retrieving known eoinpanions. Tlierefore, unlike sim- 
ilar simulations described in Kraus et al. (2005), we do 
not use DAOFIND to identify sources, but instead re- 
port the source positions directly to ALLSTAR for fit- 
ting. MHO-Tau-8 and V410-Xray3 were the only targets 
which produced significant fits for two point sources, so 
these limits should be appropriate for the detection of 
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i'-z' (mag) 

Fig. 4. — An i' vs i' — z' color magnitude diagram. Taurus 
members are shown as filled circles; objects which are classified as 
background stars based on their position in this CMD are shown 
as open circles. The SDSS field main sequence (solid line), 1 and 
2 Myr isochroncs (dashed lines), and the detection limits of the 
survey (dot-dashed line) arc also shown. The components of both 
binary systems are connected with a dotted line; the candidate 
secondary component of V410-Xray3 is well off the right side of 
the plot. The i' photometry for V410-Xray3 B is highly uncertain, 
so its i' — z' color is not likely to be accurate. The statistical error 
bars for binary components arc associated with each point; the 
error bars for single objects are smaller than the symbols, so they 
are not shown. 



new companions as well. 

In Figure 5, we show the detection limits for the four 
representative targets in V and z', as a function of sep- 
aration, at which we can detect > 10%, 50%, and 90% 
of the companions. At small separations (<20 AU), the 
50% detection thresholds roughly scale with the bright- 
ness of the primary; the survey limits are similar for all 
objects in terms of Am. The detection limits converge 
to constant values at large separations; in the case where 
the noise is background-dominated, our 5cr detection lim- 
its are z'=23.8, i' =25.1, and V=26.7. The simulations 
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demonstrate that we potentially could identify bright, 
equal-mass pairs as close as 1 pixel (0.027"; 4 AU; ^0.5 
A/D) and binaries with mass ratio g = 0.1 (Az'—A) at 
>10 AU. We also show the locations (in Am and separa- 
tion) of the companion to MHO-Tau-8 and the candidate 
companion to V410-Xray3 (for z' only). These results in- 
dicate that the probability of detecting close binary com- 
panions in the z' images is near unity for even the faintest 
systems, but the probability of detecting similar systems 
in the V filter declines rapidly for the faintest targets. 
Our inability to resolve the V410-Xray3 system in the 
V band suggests a limit oi V > 20.3 for the secondary 
component; the corresponding color limit {V — z' > 4.7) 
is consistent with the spectral type (M7.7) we determine 
in Section 4.3. 

Our simulations suggest that it is not easier to resolve 
extremely red targets at shorter wavelengths; given the 
PSF stability of HST, the improved number statistics 
at redder bands outweigh the superior diffraction-limited 
resolution at bluer bands. However, the additional color 
information was useful in confirming the identify of as- 
sociation members. Also, this is the first large set of 
homogeneous observations of young VLMOs shortward 
of 6000 angstroms. As we describe in Section 4.3, this 
allows us to characterize a previously unexplored regime 
in their spectral energy distributions. 

3.4. Uncertainties in Binary Properties 

A similar Monte Carlo routine was used to test the 
uncertainties in the measurements of MHO-Tau-8 and 
V410-Xray3. We used ADDSTAR and the average PSF 
to construct 100 simulated images, given the positions 
and brightnesses reported for the primary and the sec- 
ondary for the real images, and then used ALLSTAR to 
perform PSF-fitting photometry on these simulated im- 
ages. The standard deviation in separation (^^0.5 mas) 
reported for the simulated images of MHO-Tau-8 is con- 
sistent with that calculated from the standard deviation 
in separation between the three filters (^^2 mas), and the 
standard deviation in the flux ratio Am (0.04 magnitudes 
in each filter) is consistent with that determined from the 
magnitudes reported by ALLSTAR for each object. The 
uncertainties predicted by our simulations are generally 
lower than those observed, but this is likely due to the 
small variations in PSF width observed for each target 
(Section 2.2). We found similar results for the z'-band 
images of V410-Xray3 (0.5 mas, 0.03 magnitudes), but 
ALLSTAR only found a two-source fit for this object in 
87% of the i'-band images, and the standard deviation 
among successful fits was much higher in both separa- 
tion (^7 mas) and flux ratio (0.4 magnitudes). These 
uncertainties are dominated by the measurements of the 
secondary component, and match the uncertainties from 
the observations. 

We also conducted Monte Carlo tests to determine the 
probability of mistakenly identifying a true single star 
as a binary. We constructed a series of simulated im- 
ages (100 each for several objects spanning our sample's 
range of brightness), and then tried to fit each object 
with two point sources. We found that this never pro- 
duced a successful fit, so the probability of an erroneous 
binary identification due to statistical errors appears to 
be low (< 1%). However, since the PSF width varies 
between observations, there is also a systematic compo- 



nent to the uncertainty in binary values. We have used 
the average PSF for all fits, but since there are no other 
bright objects in either binary field, it is not possible to 
determine whether our wide or narrow PSFs might be 
more appropriate in a specific case. Experiments show 
that the separation and position angle for each binary 
are independent from the choice of template PSF, but 
the photometry is not. If we attempt to fit the binaries 
with the wider template PSFs, then ALLSTAR is unable 
to produce a significant fit. If we attempt to fit the bina- 
ries with the narrower template PSFs, ALLSTAR reports 
values for Am which are systematically lower by ~50%. 
This suggests that higher-resolution observations will be 
required to determine accurate component fluxes. 

4. ANALYSIS 
4.1. VLMO Binary Frequency 

For our 22 M5.75-M9.5 targets, we have confirmed a 
previously-identified binary VLMO system and identified 
a new candidate binary, both with separations of ^6 AU. 
The observed binary fraction is thus either 4.5±4.5% or 
9.1±6.4% for separations >4 AU. As can be seen in Fig- 
ure 5, the detection of faint companions is difficult at 
separations comparable to the PSF width (58 mas in V). 
Consequently, the total binary fraction in this separation 
range may be higher. However, studies of VLMO bina- 
ries in the field and in Upper Scorpius have found that 
companions preferentially have mass ratios q > 0.6 de- 
spite the surveys' sensitivity to companions with lower 
mass ratios (Close et al. 2003; Bouy et al. 2003; Kraus 
et al. 2005), which corresponds to Az < 1 according to 
the models of Baraffe et al. (1998). If Taurus VLMOs 
have a similar mass ratio distribution, then the detection 
thresholds in Figure 5 predict that we should identify all 
VLMO companions at separations >4 AU in at least z' 
and likely in the other filters as well. Even if these lim- 
its are optimistic, the identification of V410-Xray3 as a 
candidate binary near the minimum measurable separa- 
tion, ^5 AU, implies that the survey should be complete 
to Az' = 2 (q=0.45) for separations >4 AU. It is there- 
fore unlikely that there are many additional binaries at 
separations >4 AU with q >0.45 in this sample. 

The intrinsically higher luminosity of binaries makes 
them easier to identify in membership surveys, which can 
also bias the binary frequency to larger values. However, 
MHO-Tau-8 and V410-Xray3 are substantally brighter 
than the detection limits of their discovery surveys and 
would have been easily detected without the additional 
flux of their companions, so this effect should not bias 
our multiplicity results. 

4.2. Limits on Planetary- Mass Companions 

The high dynamic range of our deep exposures would 
allow us to directly image wide planetary-mass compan- 
ions. In Figure 5, we indicate the predicted brightness 
of some representative masses of planet ary companions 
based on the 1 Myr DUSTY models of iChabrier et alJ 
(^000). These values assume an extinction of Av~1.5, a 
typical value for objects in Taurus. The extremely red 
colors predicted for planetary-mass objects imply that 
z' observations provide the strictest limits on planetary 
companions. Based on the lack of detections, we con- 
clude that there are no planetary companions with mass 
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Fig. 5. — Five-sigma detection frequencies (10%: solid, 50%: dashed, and 90%: dash-dotted) as a function of separation in each filter 
for single VLMOs of maximum, intermediate, and minimum brightness: GM Tau, KPNO-Tau-2 and KPNO-Tau-7, and KPNO-Tau-4, 
respectively. Corresponding brightnesses of potential planetary-mass companions (assuming Ay = 1.5) are shown on the right for the z' 
plots. The brightness of the primary object is denoted with a star to allow conversion to Am values, and the vertical dotted line indicates 
a separation of 20 AU at the distance of Taurus (145 pc). The filled circles mark the separation and Am values for the binary companions 
to MHO-Tau-8 and V410-Xray3. 



>3 Mj at projected separations larger than 280 mas (40 
AU) for the brightest targets or 140 mas (20 AU) for the 
faintest targets. For comparison, these limits would have 
allowed for a significant detection of the planetary-mass 
companion to 2MASSWJ 1207334-393254, a substellar 
member of the TW Hya association, which has a pro- 
jected separation of ^^54 A U, a flux ratio of Az'^ 7, and 
a predicted mass of -5 Mj (|Ch auvm eTaOllOOl. Other 
surveys in Cham I (12 targets; Neuhauser et al. 2002) 
and IC 348 (37 targets; Luhman et al. 2005a) which in- 
cluded VLMOs in the same mass range (<0.12 Mq) also 
found no planetary-mass companions. 

4.3. Inferred (Sub) stellar Properties 
4.3.1. Masses and Spectral Types 



In Table 4, we give the inferred spectral types and 
masses for all of the VLMOs in our sample. Spectral 
types for single VLMOs and for binary primaries in our 
sample are taken from the discovery sources listed in 
Table 1 and were determined via low- or intermediate- 
resolution spectroscopy. The masses for this sample are 
estimated from the 2-Myr mass- magnitude-temperature 
relations of Baraffe et al.l lll99 8fl and the temperature- 
SpT relations of L uhman et alJ l)2003lJl . and range from 
0.015 to 0.14 A/0. Large systematic errors may be 
present in these and all pre-main sequence models (e.g. 
Baraffe et al. 2002; Hillenbrand & White 2004; Close et 
al. 2005; Reiners et al. 2005), so they are best used for 
relative comparison only. 

The models of Baraffe et al. predict that the mass ra- 
tio {q = rris/mp) and difference in component spectral 
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TABLE 4 
VLMO (Sub) Stellar Properties 



Name 


SpT'^ 


M (Mof 


Av 


L/Lms 


log L (Lq) 


R(i?0) 


Accretor?'^ 


E{V - Vms) 


CFHT-Tau-1 


M7 


0.055 


4.8 


54 


-1.4 


0.73 


N 


-0.59 


CFHT-Tau-2 


M7.5 


0.04 


2.7 


35 


-1.6 


0.56 


N 


-0.35 


CFHT-Tau-3 


M7.75 


0.035 


1.8 


29 


-1.8 


0.50 


N 


-0.71 


CFHT-Tau-4 


M7 


0.055 


4.4 


210 


-0.8 


1.45 


Y? 


-0.61 


KFNO- iau-1 


M8.5 


0.02 


1.5 


10 


-2.3 


0.27 


N 


0.08 


KPNO-Tau-2 


M7.5 


0.04 


0.7 


17 


-2.0 


0.39 


N 


-0.25 


KPNO-Tau-3 


M6 


0.1 


2.0 


18 


-1.6 


0.47 


Y 


-0.74 


KPNO-Tau-4 


M9.5 


0.01 


4.0 


31 


-2.0 


0.42 


Y? 


-1.18 


KPNO-Tau-5 


M7.5 


0.04 


0.3 


51 


-1.5 


0.68 


N 


-0.35 


KPNO-Tau-6 


M8.5 


0.02 


0.7 


9 


-2.3 


0.26 


Y 


-0.95 


KPNO-Tau-7 


M8.25 


0.017 


0.1 


11 


-2.2 


0.29 


Y 


-0.12 


KPNO-Tau-8 


M5.75 


0.11 


0.7 


15 


-1.7 


0.44 


N 


-0.02 


KPNO-Tau-9 


M8.5 


0.02 


1.7 


7 


-2.5 


0.22 


N 


0.43 


KPNO-Tau-12 


M9 


0.014 


1.6 


3 


-2.8 


0.14 


Y 


-2.15 


KPNO- iau-14 


M6 


0.1 


4.2 


123 


-0.8 


1.22 


N 


-0.42 


TV T T T nn A 

MHO-Tau-4 


M7 


0.055 


0.8 


126 


-1.1 


1.12 


N 


-0.39 


MHO-Tau-5 


M7 


0.055 


0.3 


191 


-0.9 


1.38 


Y? 


-0.39 


MHO-iau-8 AB 


Md+Md.6 


0.17 


1.7 


158 


-0.7 


1.38 


N 


-0.38 


LH 0419+15 


M7 


0.055 


1.2 


12 


-2.1 


0.35 


- 


-0.22 


V41U Xray-o AB 


TV ire^ 1 TV (Tf^ T 

M6+M7.7 


0.14 


1.5 


81 


-1.0 


0.99 


N 


-0.54 


V410 Anon- 13 


M5.75 


0.11 


5.6 


60 


-1.0 


0.87 


Y 


-0.57 


GM Tau 


M6.5 


0.08 


0.2 


49 


-1.1 


0.74 


Y 


-1.23 


USco 55 


M5.5 


0.010 


0.1 


18 


-1.5 


0.49 


N 


0.18 


USco 66 


M6 


0.07 


-0.4 


16 


-1.7 


0.44 


N 


-0.58 


USco 67 


M5.5 


0.10 


0.4 


18 


-1.5 


0.49 


N 


0.25 


USco 75 


M6 


0.07 


0.2 


19 


-1.6 


0.48 


N 


-0.43 


USco 100 


M7 


0.05 


-0.2 


34 


-1.6 


0.58 


N 


-0.40 


USco 109 


M6 


0.057 


0.4 


9.1 


-1.9 


0.33 


N 


-0.33 


USco 112 


M5.5 


0.10 


-0.2 


6.8 


-1.9 


0.30 


N 


-0.63 


USco 128 


M7 


0.05 


0.5 


10 


-2.2 


0.32 


N 


-0.28 


USco 130 


M7.5 


0.04 


0.4 


12 


-2.1 


0.33 


N 


-0.28 


USco 131 


M6.5 


0.06 


1.3 


6.7 


-2.2 


0.27 


N 


-0.08 


USco 132 


M7 


0.05 


1.9 


14 


-2.0 


0.37 




-0.40 


USco 137 


M7 


0.05 


1.0 


3.3 


-2.6 


0.18 




-0.28 



^Spectral types for binary secondaries are inferred from photometry presented here; others are from their discovery source (Section 4.3) 
'^Masses are determined from the models of Baraffe et al. (1998). 
•^Muzerolle et al. 2003; Mohanty et al. 2005; Muzerolle et al. 2005. 



types (ASpT) for young VLMO binaries are a function 
of the primary-to-secondary flux ratios Am with only a 
minor mass dependence; we report these quantities as 
determined from the flux ratio Az' in Table 3 for the 
detected binary VLMOs, MHO-Tau-8 and V410-Xray3. 
The models do not report z' magnitudes so we calcu- 
lated these from the model Ic magnitudes, the Ic — i' 
transformations determined in Section 2, and the %' — z' 
colors for the appropriate spectral type. The uncertain- 
ties reflect the uncertainties in the flux ratios and the 
scatter in this relation in the models, but do not include 
any systematic uncertainties from the models. The val- 
ues determined from AV and Ai' for MHO-Tau-8 are 
roughly consistent with these results, but the values for 
Az' are the only reliable result for V410-Xray3, so we 
use Az' in all cases for uniformity. 

4.3.2. Luminosities and Extinctions 

Our results, in combination with other surveys, also 
allow us to determine the luminosity and extinction for 
each VLMO by fitting the obsca'vcd spectral energy dis- 
tribution (SED) with a reddened dwarf SED of the same 
spectral type. This allows us to test for potential opti- 
cal excesses that previous observations at longer wave- 
lengths would not have detected. The VLMO SEDs are 
constructed from Vi'z'JHK photometry from this work 



and from the Two Micron All-Sky Survey (2MASS; Cutri 
et al. 2003). Our standard SEDs are constructed from 
the mean i'z'JHK photometry for field stars in the SDSS 
(West ct al. 2005) and 2MASS (Lcggett et al. 2002), plus 
the mean V — K colors reported for the field by Reid et al. 
(2003). Typical uncertainties are ;$0.05 magnitudes for 
our observations and ~0. 1-0.2 magnitudes for the field 
star photometry. Any significant difference between a 
VLMO SED and the corresponding field SED should be 
a result of the intrinsically higher luminosity of our tar- 
gets (pre-main sequence VLMOs have larger radii) or 
extinction due to obscuring material. These will result 
in a constant multiplication and a wavelength-dependent 
multiplication to the fiux in each band, respectively. We 
solve for these constants, which correspond to L/Lms 
and Av, by a least-squares fit of the set of equations: 

miyLMO = rUiJield + T^^V " '^■^lOQ- 

Av LmS 

where i denotes five filters used in the fit (i'z'JHK) 
and the reddening coefficient for each band is taken from 
Schlegel et al. (1998). Variations in distance will also 
be included in the luminosity term, but if the depth is 
similar to the apparent width of Taurus- Auriga on the 
sky (~10 deg; 25 pc), these should be relatively minor. 

We list the inferred luminosity ratio, total luminosity, 
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radius, and extinction for each Taurus VLMO in Table 
4. The total luminosity is determined from the K-band 
magnitude, corrected for extinction, and the bolometric 
corrections of Leggett et al. (2002), which are BCx~3.1 
for mid to late M dwarfs. The radius is determined 
from the luminosity ratio {L/Lms k {R/ RmsY) and 
the main sequence mass-radius relations of Baraffe et 
al. (1998), which find that Rms/Rq Mms/Mq for 
M5-M9 dwarfs at ages >1 Gyr. We also performed some 
experiments in varying the spectral type of the field stan- 
dard; these imply that the uncertainty in spectral type 
(typically half a subclass) corresponds to uncertainties 
of ~0.5 in Av, ~20% in L/Lms, and ~10% in R. The 
uncertainty in the bolometric correction is typically 
for spectral types later than M5 (Leggett et al. 2002). 
The statistical uncertainties from the fitting process are 
generally insignificant compared to these systematic un- 
certainties. These experiments also reveal a degeneracy 
in the determination of extinction and luminosity; fitting 
with an earlier or later spectral type yields systematically 
different extinctions and luminosities, but the goodness 
of fit does not decrease significantly for small changes (<1 
subclass). The same fitting process has been applied to 
the binaries by considering the total system brightness 
in each filter and assigning it the spectral type of the 
primary fitting it with a field standard corresponding to 
the spectral type of the primary. However, since the ob- 
served spectrum is a composite of both objects' spectra, 
the results may be biased. 

In Figure 6, we present SEDs for the 18 Taurus VLMOs 
which produced successful fits to the dwarf SEDs (within 
2(t). Four stars could not be fit well and are discussed 
below. The solid lines and filled circles denote VLMO 
SEDs which have been corrected for extinction and the 
dashed lines denote the best fit dwarf SEDs which have 
been shifted upward by log{L/ Lms) to correct for their 
intrinsically lower luminosity. Typical uncertainties in 
the SED flux measurements are ~0.02 dex for the VLMO 
SEDs, -0.05 dex for the r'i'z'JHK fluxes in the dwarf 
SEDs, and ~0.15 dex for the V band fluxes in the dwarf 
SEDs. Since our results typically agree with the dwarf 
SEDs to within the uncertainties, the common assump- 
tion of dwarf colors for pre-main sequence objects ap- 
pears to be valid for this age and mass regime. However, 
as we discuss in Section 4.3.3, the targets appear to be 
systematically brighter than equivalent dwarf SEDs in 
V, yielding bluer optical colors. In Figure 7, we present 
SEDs for the four VLMOs with signiflcant discrepancies 
in their flt: KPNO-4, KPNO-6, KPNO-12, and GM Tau. 
The first three objects deviate significantly only in their 
V band measurement, but the optical and near-infrared 
portions of the SED for GM Tau appear to be inconsis- 
tent; we can not produce a flt to its dwarf reference by 
any single combination of Ay and L/Lms- Based on the 
slope of its NIR SED, GM Tau appears to have been sig- 
niflcantly fainter and/or possessed a significant K-band 
excess during its epoch of observation in 2MASS. The ex- 
cellent agreement between the optical and near-infrared 
values for all other objects, which were taken at different 
epochs, suggest that they are typically not variable at 
amplitudes of >0.2 magnitudes. 

4.3.3. Optical Excesses 



Nine of the VLMOs in our sample have been observed 
to possess spectroscopic signatures of active accretion of 
circumstellar material by Muzerolle et al. (2003), Mo- 
hanty et al. (2005), and Muzerolle et al. (2005). Ac- 
creted material falling onto more massive T Tauri stars 
is typically heated to higher temperatures than the stel- 
lar photosphere, leading to an optical excess (Basri & 
Batalha 1990; Hartigan et al. 1991). However, since the 
rate of accretion appears to depend on system mass (e.g. 
Muzerolle et al. 2003), the low mass accretion rates of T 
Tauri VLMOs should not produce any significant excess. 
Nonetheless, as we note in Section 4.3.2, several VLMOs 
show signs of an optical excess which might be a result 
of accretion. In Table 4, we identify the known accretors 
and non-accretors and list the V-band excess or deficit 
(in magnitudes) observed for each target. The four ob- 
jects with the most signiflcant V band discrepancies are 
conflrmed or probable accretors. 

In Figure 8, we show a plot of the V-band excesses 
and deficits as a function of spectral type for each tar- 
get. Filled circles denote Taurus VLMOs; symbols for 
known accretors (based on Ha or other line profile diag- 
nostics) are larger than those of nonaccretors. Open cir- 
cles denote Upper Sco VLMOs, none of which are known 
accretors. The dotted line denotes the V-band excess 
one would mistakenly infer if the color were infiuenced 
by surface gravity effects, estimated here using the in- 
trinsically bluer V — I colors of M5-M8 giants (Bessell & 
Brett 1988). Most of the known accretors sit preferen- 
tially higher in the plot relative to the non-accretors and 
relative to the giant locus, implying that they possess 
intrinsic excesses which are related to accretion or to the 
assessment of their accretion status. 

Even non-accretors seem to possess small V-band ex- 
cesses over typical dwarf colors, though they are not as 
high as typical giant excesses and do not depend signifi- 
cantly on spectral type; this suggests that young VLMOs 
tend to be bluer than dwarfs of the same spectral type. 
The mean and standard deviation for this discrepancy 
among known non-accretors are 0.28±0.33 magnitudes; 
in Figure 8, the mean value is denoted with a solid line 
and the +2a limit with a dashed line. If we regard this 
mean value as a measure of the true photospheric flux 
for young VLMOs, then the four objects identifled ear- 
lier still possess excesses signiflcant at >+2a. 

In Table 5, we present the (sub)stellar properties for 
our Upper Sco sample (Kraus et al. 2005), which we de- 
rive using the methods described in Section 4.3. There 
are no known accretors, and the mean V-band excess 
(0.27±0.08) is consistent with the mean excess for nonac- 
cretors in Taurus (0.28±0.10). As we show in Figure 8, 
all Upper Sco VLMOs fall below the +2a criterion for a 
statistically signiflcant optical excess over the nonaccre- 
tor value. This supports the suggestion that accretors 
may show optical excesses over nonaccretors, though the 
older age of Upper Sco may also play a role. 

5. DISCUSSION 

5.1. (Suh)stellar Multiplicity 

Multiplicity surveys of VLMOs have suggested a 
fairly uniform set of binary proper t ies in the fleld 
JClosc ct al." 2003; Bo uv et all I2003t IBurgasse r et al. 
[2003^: ..Siegler et al. 2005i) . in the Pleiades (,,Martin et al. 
12003ft . and in the young OB Association Upper Scorpius 
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Fig. 6. — Spectral energy distributions for the VLMOs in our Taurus sample which produced successful fits. Solid lines and filled circles 
denote the target SEDs which have been corrected for extinction. Dashed lines denote field SEDs which have been corrected for the higher 
luminosity of our targets. 



(Kraus et al. 2005). No companions with wide separa- 
tions (>20 AU) or with unequal mass ratios {q < 0.6) 
were found, despite sufficient sensitivity for their de- 
tection. Several wider systems have been discovered 
serendipitously or during ongoing surveys (Gizis 2000; 
Luhman 2004a, 2005c; Phan-Bao et al. 2005; Billeres et 
al. 2005). However, the corresponding survey statistics 
are not available, so it is difficult to determine a fre- 
quency for wide VLMO binaries. The binary fractions 
observed were also significantly l ower than the binary 
fracti ons of 57% for field G-dwarfs ('Duau ennov fc Mavorl 
1991 and 35-43% for field M-dwarfs (R.ci d fc Gizislll997t 
Fischer fc Ma^\im^ . We summarize these results in 
Table 5. 

The study presented here is the first large survey of 
young VLMOs in Taurus-Auriga with sufficient resolu- 
tion to identify binaries at separations seen for field VL- 
MOs. Our confirmed binary fraction, 4.5 ±4.5% for sep- 



TABLE 5 

Results for Previous Multiplicity Surveys 



Survey 


SpTs 


Minimum 


Binary 






Sep (AU) 


Frequency 


Close et al. (2003) 


M8.0-L0.5 


3 


15±7% 


Bouy et al. (2003) 


M8.0-L8.0 


1 


18.8±3.7% 


Burgasser et al. (2003) 


T5.0-T8.0 


1 


9+1% 


Siegler et al. (2003) 


M6.0-M7.5 


3 


9tt% 


Martin et al. (2003) 


M6.0-M9.5'' 


7 




Kraus et al. (2005) 


M5.5-M7.5^ 


5 


25+i'5% 


Current Survey 


M5.75-M9.5'' 


5 


9+l°% 



^Thc surveys by Martin et al. and Kraus et al. observed BDs in the 
Pleiades and Upper Sco, respectively. Since these regions are young, 
the spectral type range corresponds to lower masses than in the field. 
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Fig. 7. — Spectral energy distributions for the four VLMOs in our sample which produced anomalous fits. 
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Fig. 8. — V-band excesses as a function of spectral type. Filled circles denote Taurus members, open circles denote Upper Sco members, 
and large symbols denote confirmed or probable accretors. A solid line shows the mean value for all non-accretors in our sample, and the 
dashed line shows the corresponding 2(7 limit for the existence of significant optical excesses. The four objects with significant excesses are 
labeled. 
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arations >4 AU, is lower than, but marginally consistent 
with the results cited for the field. If the binary candi- 
date V410-Xray3 is confirmed, the binary fraction in this 
separation range would then be 9.1 ± 6.2%, which falls 
within the range of results seen in the field. However, 
the field surveys are sensitive to companions with sepa- 
rations as small as ~ 1 AU and found many companions 
within 4 AU that we could not have detected because of 
the larger distance to our targets. Our binary fraction 
for companions with a > 4 AU only sets a lower limit 
on the total binary fraction for VLMOs. The mass ratio 
and separation for MHO-Tau-8 (q=0.7, a~6 AU) and the 
separation for V410-Xray3 (~6 AU) are consistent with 
the field, but the mass ratio for V410-Xray3 may be sub- 
stantially less. However, this result is based only on the 
z' observations and should be confirmed with foUowup 
observations. 

The methodology of this survey is identical to our Up- 
per Scorpius survey and their similar distances imply 
similar spatial resolution, so the two results can be com- 
bined to produce a total VLMO binary fraction (either 
4 or 5 of 34, corresponding to 12 ± 6% or 15 ± 7%) which 
is also consistent with the field. Our binary fraction for 
Taurus (9±5%) is somewhat lower than that for Upper 
Scorpius (25±16%), though the difference is not statisti- 
cally significant (<2ct). However, larger surveys (partic- 
ularly in Upper Sco) will be required to test whether the 
binary fraction depends on mass. 

5.2. The Mass Dependence of Multiplicity 

The distinct binary properties of field BDs, relative 
to those of stars, indicate that these properties are 
mass dependent. H owever, the f o rm o f this depen- 
dence is not known; IKroupa et al.l lj2003(l . I Close et "all 
l)2003j) . and Kraus et al. (2005) interpret current data as 
a sharp transition near the stellar/substellar boundary 
while "tuhman <j2004bj) and White et al. (2006) argue 
for a smooth mass dependence. The field binary fraction 
appears to decrease from 57% for G dwarfs to 35-43% for 
M dwarfs (Duquennoy & Mayor 1991; Fischer & Marcy 
1992; Reid & Gizis 1997), but the trend with mass is 
only marginally significant. 

Our combined young VLMO sample for Taurus and 
Upper Scorpius is large enough to investigate a possible 
mass dependence in the VLMO binary frequency over 
the single order of magnitude mass range of our survey 
(0.015-0.12 Mq). All five binaries and candidate binaries 
were found among the 13 most massive VLMO targets 
{Al > 0.07 Mq); none were found among the 21 lower- 
mass VLMOs. If the probability of binarity is 5/34 for 
all members of our sample, then the probability that all 
five binaries will fall among the 13 most massive tar- 
gets is (13!/8!)/(34!/29!) - 0.005. The possible double- 
lined spectroscopic binary KPNO-Tau-14 (Mohanty et 
al. 2005) also falls in the upper mass range, though its 
separation is well inside our survey limit. The detec- 
tion limits in Figure 5 demonstrate that bright binaries 
are more easily identified due to superior photon statis- 
tics, but at least the three high-confidence binaries in the 
combined sample (USco-55, USco-66, and MHO-Tau-8) 
could have been identified around even the faintest tar- 
gets. Also, since this trend is seen at lower significance in 
both regions, it should be robust against variations in ini- 
tial conditions. The implication is that the binary frac- 



tion may decline over the mass range, though a constant 
binary fraction combined with a shift toward smaller sep- 
arations could produce the same effect in our resolution- 
limited sample. 

We can further investigate the mass dependence of 
multiplicity in Taurus by combining our results with 
a binary census of more massive Taurus members by 
White et al. (2006). This speckle interferometry sur- 
vey included 44 members of Taurus in the mass range 
0.06 < M < 1.5 Mq and, in combination with previ- 
ous results for 41 additional targets (Leinert et al. 1993; 
Ghez et al. 1993; Simon et al. 1995; Duchene 1999), 
identified 23 binaries with separations > 27 mas (4 AU). 
The sample is reported to be complete for separations 
of 9-460 AU and mass ratios >0.1 (Ai^<2.75), which 
falls within the completeness range of our survey (Sec- 
tion 3.3). They report that the binary frequency ap- 
pears to peak for stars with masses of 0.3-0.7 Mq and 
decline toward lower masses, though the results are not 
statistically significant (<2cr). Our VLMO sample sup- 
ports this mass dependence, since we found no binaries 
with separations in the completeness regime. However, 
as discussed in Kraus et al. (2005), a similar survey in 
Upper Scorpius by Kohler et al. (2000) found no evi- 
dence of a similar decline until immediately above the 
stellar/substellar boundary (0.1-0.2 Mq). This implies 
that the mass dependence of multiplicity could also have 
a regional dependence. 

5.3. A Future Dynamical Mass for MHO-Tau-8? 

Theoretical models of the mass-luminosity- 
temperature relation are still largely uncalibrated 
for low masses and young ages (e.g. Baraffe et al. 2002). 
Surface gravit y measurements for young brown dwarfs 
in Upper Sco ijMohantv et al.ll2'00,5D have provided some 
evidence that the theoretical models underestimate 
masses for very low-mass stars and brown dwarfs and 
overestimate masses for the least massive brown dwarfs, 
but the measurement of dynamical masses for low-mass 
binaries will be required to directly calibrate the models. 
To date, this has been done only for three young 
VLMOs, the older (-50-100 Myr) brown dwarf AB Dor 
C (Close et al. 2005; Luhman et al. 2005b) and both 
components of the young substellar eclipsing binary 
2MASS J05352184-0546085 (Stassun et al. 2006). 

Relative to the speckle interferometry measurement by 
White et al. (2006), MHO-Tau-8 B has traced -40 de- 
grees of its orbit around MHO-Tau-8 A within only —5 
years; this suggests an orbital period of —45 years. If 
we assume a circular face-on orbit, then we estimate 
from Kepler's Law and the system parameters [a > 6 
AU, P — 45 years) that the total mass of the system 
is >0.11 Mq, which is consistent with the total system 
mass implied by the models of Baraffe et al. (1998), 
—0.17 Mq. Since this lower limit is substantially lower 
than the predicted total mass, the orbit of the secondary 
is likely either inclined relative to the plane of the sky 
(implying a larger semi-major axis) or highly eccentric 
and near apastron. In either case, additional astrometric 
measurements for this system should yield precise orbital 
parameters and better limits on the dynamical masses of 
its components over the next decade. 

5.4. VLMO Optical Excesses 
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The rate of (sub)stellar mass accretion onto T Tauri 
stars is typically estimated either by measuring the opti- 
cal continuum excess or by modeling its effect on emission 
line features that result from outflows as accreting ma- 
terial impacts the photosphere. Previous observations 
of brown dwarfs and low- and intermediate-mass stars 
have found mass accretion rates that scale roughly as 
M oc M2 for masses of 0.03-3.0 Mq (Gulbring et al. 
1998; White & Ghez 2001; Muzerolle et al. 2003; Calvet 
et al. 2004). Specifically, these surveys found mass accre- 
tion rates ranging from '^lO^^ MQ/yr for intermediate- 
mass T Tauri stars with mass 3 Mq to ^10^"'^"'^ M0/yr 
for T Tauri-age brown dwarfs. Recent results from Mo- 
hanty et al. (2005) and Muzerolle et al. (2005), who 
measure accretion rates by modeling of Ca II and Ha 
emission respectively, support this trend. They foimd 
that many low-mass brown dwarfs (including KPNO-6, 
KPNO-12, and possibly KPNO-4) have apparent mass 
accretion rates of ~5xl0~^^ Mq/jt. 

These low accretion rates woiild result in negligible op- 
tical excesses {E{V — Vms) < 0.01; e.g. Muzerolle et al. 
2003); the very large excesses observed for our targets im- 
ply that either these objects have much higher accretion 
rates or the V-band excesses are not the result of accre- 
tion. The mass accretion rates that would be required for 
each known accretor to generate the optical excess can 
be calculated using the method described in GuUbring et 
al. (1998) whereby the excess V-band flux is converted 
to a total accretion luminosity. This luminosity is then 
assumed to result from mass infall from the inner edge of 
a circumstellar disk onto the surface of the star via mag- 
nctosphcric accretion. Specifically, we assume that the 
radius of the inner disk edge is 5Rvlmo and the entire 
change in potential energy is radiated away via the op- 
tical excess. To correct our V-band excess luminosity to 
the total accretion luminosity, we adopt a bolometric cor- 
rection factor of 10, which is 10% less than the bolometric 
correction adopted by White & Hillenbrand (2004) for a 
slightly redder, but narrower band-pass (0.60 - 0.65 /xm). 
As emphasized by White & Hillenbrand, these bolomet- 
ric corrections are still highly uncertain, but by adopt- 
ing a value consistent with previous assumptions, we can 
minimize systematic errors in comparisons with previous 
calculations. Unfortunately, there are no photometric or 
spectroscopic data available at wavelengths shorter than 
6000 angstroms for young BDs, so our results can not be 
directly compared to any in the literature. 

The resulting mass accretion rates for the 4 stars with 
significant excesses (10^^ to 10^^ MQ/yr) arc 3-4 or- 
ders of magnitude higher than the values found via other 
methods or inferred from the M-M relations observed for 
higher-mass T Tauri stars. Given this discrepancy, we 
consider the possibility that the excess results from our 
analysis methods or from other physical processes. How- 
ever, it is statistically unlikely that 4/22 objects would 
possess >+2a excesses, or that KPNO-12 would possess 
a 6(J excess, if these accretors are part of a random dis- 
tribution with the same mean and standard deviation as 
we observe for nonaccretors. Since 4/9 accretors possess 
excesses significant at >+2a and 8/9 possess excesses 
significant at >+1(t, it appears that the V-band excesses 
are related to accretion diagnostics. Since three of the 
outliers fall among the 5 latest objects, there could also 



be a gravity-based V-band excess inherent to all late- type 
young VLMOs. However, the corresponding trend for gi- 
ants (which possess even lower surface gravity) appears 
to decline for spectral types later than M7 (Figure 8), 
and there is no sign of a similar excess for the two non- 
accreting late-type objects, KPNO-1 and KPNO-9. Fi- 
nally, the possibility of erroneous spectral classifications 
cannot explain this result because there are significant 
discrepancies in the SED fits for errors of more than one 
subclass, and the degeneracy between spectral type, Ay, 
and L/Lms that we note in Section 4.3.2 acts to preserve 
most of the inferred V-band excess under the assumption 
of a different spectral type. 

One explanation that we cannot rule out is the pres- 
ence of chromospheric activity or flaring. Stellar flares 
are characterized by transient regions with higher tem- 
peratures than the surroimding photosphere (<;10"'' K); 
they cause a temporary increase in the total luminosity, 
with much higher amplitude of variability at short wave- 
lengths. Surveys for x-ray emission from young BDs (e.g. 
Preibisch et al. 2005) have found that 1/3 to 1/2 of aU 
objects may possess significant activity. These regions 
also typically produce Ha emission; the width of the base 
of the Ha emission line is a typical test for accretion, so it 
is possible that a large flare or significant chromospheric 
activity could broaden the line base, mimicking the sig- 
nature of accretion. 

6. CONCLUSIONS 

We present the results of a high-resolution imaging sur- 
vey of 22 brown dwarfs and very low mass stars in the 

nearby (^^145 pc) yoimg (^^1-2 Myr) low-density star- 
forming region Taurus- Auriga. This survey confirmed 
the binarity of MHO-Tau-8 and discovered a new can- 
didate binary system, V410-Xray3, resulting in a binary 
fraction of 9 ± 5% at separations >4 AU. Both binary 
systems arc tight (<10 AU) and they possess mass ra- 
tios of 0.75 and 0.46, respectively. The binary frequency 
and separations are consistent with very-low-mass binary 
properties in the field, but the mass ratio of V410-Xray3 
is among the lowest known. The binary frequency and 
separations are also consistent with the trends in multi- 
plicity statistics observed for higher-mass Taurus mem- 
bers, which suggest a gradual decline in both properties 
toward low masses; the implications are that the dis- 
tinct binary statistics of low- mass systems are set in the 
formation process and that this formation process is sim- 
ilar to that which creates low-mass stars. These objects 
and a similar set in Upper Scorpius reveal another possi- 
ble mass dependence of multiplicity; all five binaries and 
candidate binaries fall in the more massive half of the 
combined sample, implying either a decline in frequency 
or a shift to smaller separations for the lowest-mass bi- 
naries. We also combine the survey detection limits with 
the models of Chabrier et al. (2000) to show that there 
are no planets or very low-mass brown dwarfs with sep- 
arations > 3M,/ at projected separations >40 AU orbit- 
ing any of the Taurus members in our sample, implying 
that planetary-mass companions at wide separations like 
2M1207b are found with frequency <5%. Finally, we ob- 
serve significant optical excesses in the spectral energy 
distributions of most targets and conclude that the tar- 
gets with spectroscopic signatures of accretion possess 
larger optical excesses than other young brown dwarfs 
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and very low mass stars. 
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ful critique of this paper. This work is based on obser- 
vations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 
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cessing and Analysis Center/California Institute of Tech- 
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